CMB Lensing Constraints on Dark Energy and Modified Gravity Scenarios 
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Weak gravitational lensing leaves a characteristic imprint on the cosmic microwave background 
temperature and polarization angular power spectra. Here we investigate the possible constraints 
on the integrated lensing potential from future CMB angular spectra measurements expected from 
Planck and EPIC. We find that Planck and EPIC will constrain the amplitude of the integrated 
projected potential responsible for lensing at 6% and 1% level, respectively with very little sensi- 
tivity to the shape of the lensing potential. We discuss the implications of such a measurement in 
constraining dark energy and modified gravity scalar-tensor theories. We then discuss the impact of 
a wrong assumption on the weak lensing potential amplitude on cosmological parameter inference. 

PACS numbers: 98.80.-k 95.85.Sz, 98.70.Vc, 98.80.Cq 



I. INTRODUCTION 

Due to extremely sensitive observations with satellite 
(0), ground-based [Q) and balloon-borne experi- 
ments (0), the Cosmic Microwave Background (CMB 
hereafter) temperature and polarization anisotropics are 
now measured with astonishing precision. Moreover, 
these measurements are in nearly perfect agreement with 
theoretical predictions of an adiabatic cold dark matter 
(CDM) model of structure formation and most of the 
CDM cosmological parameters are now constrained at 
better than a few percent uncertainty. 

Future data expected from Planck ([^) and ground- 
based experiments such as SPT (Q) will provide even 
better measurements with cosmic variance limited maps 
of CMB temperature anisotropics down to few arcmin- 
utes angular resolution. Moreover, balloon-borne experi- 
ments such as SPIDER (0) and EBEX (|]), and possi- 
ble future satellite mission EPIC (|^) will provide further 
constraints on the polarization signal, again at the cos- 
mic variance limit of polarization measurements down to 
ten arcminute angular scale or better. 

With all this high quality data expected in the very 
near future it is definitely timely to investigate if the 
theoretical predictions have reached a similar accuracy. 

The calculation of CMB anisotropics through current 
Boltzmann solvers as CMBFAST Q or CAMB is 
made in the linear regime, where the evolution equations 
can in principle be solved to an arbitrary precision. Cur- 
rent codes show deviations at ~ 0.1% level, well below 
the future experimental accuracy (see e.g. p2[). 



However, additional uncertainties are present. The 
thermal history of primordial recombination, for exam- 
ple, is still not sufficiently theoretically determined to 
match future observations (see e.g. |p^). 

While the needed accuracy in this case relies on stan- 
dard physics and it is expected to be achievable, a more 
worrying aspect concerns the assumptions on the un- 
known dark energy component. The choice of the dark 
energy framework, indeed, clearly modifies the CMB 
anisotropy angular spectrum. It is well known that the 
dark energy density and its equation of state w defined 
as the ratio between pressure and energy of the dark en- 
ergy fluid both change the angular diameter distance of 
the sound horizon at last scattering and shift the peaks 
positions in the CMB spectra (see e.g. These pa- 

rameters can be accurately measured by the use of several 
complementary observables such as luminosity distances 
of type la supernovae ([p^). 

It has however recently become clear that models of 
dark energy or modified gravity can be conceived lead- 
ing to the same background expansion of A — CDM (and 
therefore preserving the angular diameter distance at red- 
shift z ^ 1100) but with different evolution for the per- 
turbations in the energy components (see e.g |l6), |T7| ] 
and references therein). 

The CMB angular spectra in those models may have 
small but measurable differences respect to those ex- 
pected in the A — CDM model. First, on large an- 
gular scales, since the growth of perturbations is dif- 
ferent, gravitational potentials will change differently 
with time, affecting the CMB through the so-called Inte- 
grated Sachs- Wolfe (ISW) effect |18 . Recent works (see 
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e.g. 1 19 ; have focused on this large angular scale 

anisotropy signal in order to extract information on dark 
energy or modified gravity. It is clear that one could in- 
vert the argument and ask to what extent our poor the- 
oretical knowledge of the ISW signal could affect cosmic 
parameters. Fortunately this signal is dominated by the 
primordial anisotropy and present only on the very few 
first multipolcs. In few words, the impact on parameter 
extraction of a wrong theoretical prediction for the ISW 
signal is negligible. 

However, dark energy also changes the small angular 
scale anisotropy through lensing. While the physics of 
CMB lensing is well understood, the amplitude of the 
signal is indeed connected to the growth of matter den- 
sity fluctuations with redshift. Since the growth of struc- 
ture strongly depends on the dark energy component and 
since a clear physical understanding of this component is 
lacking, it is definitely possible that the lensing amplitude 
will be different from the one expected in the cosmologi- 
cal constant scenario. 

Several recent papers have shown that the Icnscd CMB 
signal may be used to study dark energy ^ ^ . 
Here we follow a complementary approach, discussing the 
impact of an unaccounted variation in the lensing sig- 
nal, motivated either by a different dark energy model or 
modified gravity, on cosmological parameter inference. 
As we show in the next section, modified gravity, for 
example, could easily increase the rms lensing signal by 
^ 20% while leaving the primary anisotropics unaffected. 
Smaller but still sizable variations in the amplitude could 
be induced by perturbations in the dark energy fluid. We 
will show that if future data will be analyzed without 
considering the possibility of a non-standard lensing sig- 
nal, this may drastically bias the conclusions on many 
cosmological parameters. 

Our paper is organized as follows: in the next Section 
we show how dark energy perturbations and/or modi- 
fied gravity could change the amplitude of the weak lens- 
ing signal. In Section III we present our data analysis 
method and the datascts considered. In Section IV we 
present the results and in Section V we conclude with 
our conclusions. 



II. THE AMPLITUDE OF THE WEAK 
LENSING POTENTIAL IN NON STANDARD 
MODELS OF STRUCTURE FORMATION 

The lensing deflection of CMB photons depends upon 
gradients in the total gravitational potential $ -I- trans- 
verse to the line of sight to the last scattering surface , 
where $ and are defined by the perturbed Robertson- 
Walker line- element 

ds^ = [- (1 + 2^-) dr^ + (1 - 2$) dx^] , (1) 

using the notation and convention of Ref . js^ . 

The evolution of the gravitational potential can be 
expressed by a a transfer function Tq,{k,T), whereby 



^'(fc, r) = T<s,{k,T)R{k) and where R{k) is the primor- 
dial curvature perturbation. The power spectrum of the 
lensing potential is given by 
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dxSq,{k]To - x)ji{kx) 



(2) 

Here PR{k) is the primordial power spectrum, tq ~ x 
is the conformal time at which a given photon was at 
the position x?^, and the lensing source, in the standard 
scenario, is given by (see e.g. pGf): 



S^{k,TQ - x) = 2T*(fc,ro - x)jiikx) 



X* 



X 
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(3) 



The effects of lensing on temperature and polarization 
anisotropy have been extensively presented in the liter- 
ature. In the case of temperature, lensing modifies the 
damping tail, smearing the acoustic oscillations. A simi- 
lar effect but more pronounced is present in the EE po- 
larization and TE cross temperature-polarization spec- 
tra, while lensing also introduces an extra B-mode po- 
larization signal. The lensed B-modc peaks at tens of 
arcminute angular scales and its amplitude is directly 
proportional to the lensing power spectrum. It is clear 
from the expression of the lensing angular spectrum in 
Equation |^ that the lensing amplitude will depend on 
the growth of dark matter perturbations. The evolution 
of the dark matter perturbations can be altered by the 
dark energy component, affecting the weak lensing signal. 

While the equation of state of dark energy w can be 
constrained by complementary observables, the dark en- 
ergy perturbation sound speed, (see e.g. that 
probes the nature of dark energy fluctuations, is more 
elusive. It is useful therefore to investigate the impact of 
a different choice for on the weak lensing signal. Fix- 
ing w = —0.8, that is already at the border of current 
measurements that are pointing towards w = — 1±0.1, a 
^5% difference is present between the lensing potential 
with = 1 and cj. = 0. Perturbations in the dark energy 
will be present only if the dark energy equation of state 
will be different from w = — 1 and closer the equation of 
state will be to this value, smaller will be the effect. One 
could therefore argue that if future measurements will ac- 
curately determine that the equation of state is extremely 
close to —1 then there will be no expected modification 
to the weak lensing signal^. 

However, if modified gravity is at works, it is possi- 
ble to obtain a background expansion that mimics ex- 
actly A — CDM but still with a different lensing ampli- 
tude. Several modified gravity models have been pro- 
posed. Here we follow the approach presented in 
where a modified gravity changes the relations between 



^ This is not completely true since if dark energy interacts with 
dark matter then perturbations will be possible even if the, ef- 
fective, measured value of w will be close to —1 
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the two Newtonian potentials and between the same po- 
tentials and matter perturbations. Modified gravity can 
indeed be parameterized with two functions 7 and fi (see 
Hi) such that: 



2M2 



pi(a, k)pA , 



and 



(4) 



(5) 



where pA is the co-moving density perturbation and Mp 
is the Planck mass. 

The functions /x and 7 encode the effects due to mod- 
ified gravity or dark energy and can be expressed as 



and 



7(0, fc) 



1 + Xja'^ 



1 -t- A^a^ 



(6) 



(7) 



respectively, where Ai is the Compton wavelength of the 
modified force, /3i and (32 are linked to the coupling of the 
scalar field with matter and s sets the time-dependence 
of /X and 7. 

Referring to scalar-tensor theories, these parameters 
satisfy the equations 



Pi 



A? 



2- A 



A| 



(8) 



If we choose a particular value for (32 and s we 
specify a particular class of theories; for example, we 
can consider (3i ~ 4/3, (32 = 1/2 and s = 4 and work 
with a general f{R) theory |34| where Ai can be related 
to the expression of f{R) thanks to the definition of the 
mass related to the modified force = (1 + Jr)/ Jrm 



A^- ^ 



3/ 



RR 



l + fR 



(9) 



where Jr = df/dR, Jrr = d'^f/dR^ and the subscript 
denotes quantities evaluated at the present time. 

Using a modified version of the CAMB code (MG- 
CAMB, see |3^) we have therefore computed the angular 
spectra for different values of X\. 

In Figure Q we plot the Ci temperature anisotropy 
spectrum for different values of Ai, including Ai = re- 
producing the ACDM model; we also plot the relative 
difference between modified Ci and the ACDM model. 
It is possible to notice that while differences arise for 
large scale due to the ISW effect also differences are well 
present at small angular scales due to the modification 
in the lensing potential. 
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FIG. 1: CMB anisotropy power spectra (Top Left Panel) in 
a modified gravity model that mimics the background expan- 
sion of A-CDM. As we can see (Top Right Panel) the main 
differences arise on large angular scales due to the Integrated 
Sachs Wolfe effect and on small angular scales due to lensing. 
If lensing is not considered the spectra are identical on small 
scales (Bottom Panel). 



In Figure left panel, we plot the corresponding lens- 
ing potential spectra. As we can see, modified gravity 
could enhance the lensing signal in a substantial way, 
and with greater effects respect to c^. One should also 
bear in mind that the models considered mimic a cosmo- 
logical constant for the background evolution. CMB lens- 
ing could therefore be a powerful method to disentangle 
modified gravity from a cosmological constant. Values of 
\\ > 10^ Mpc^ could be ruled out by solar system tests 
even for f{R) theories that shows a chameleon mecha- 
nism; for example we can relate Ai to the parameter Jrq 
of 6^ obtaining 



which leads, for the chosen value of Af , to a |/ijo| ~ 0.1 
which is barely compatible with solar system bounds. 

Also in Figure H, right panel, we plot the lensing poten- 
tials for a Chamelon scalar-tensor theory (see [0). This 
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Chameleon 





FIG. 2: Lensing potential for diflterent f{R) models (left 
panel) and Chameleon models (right panel) in function of 
different choices of the compton length Ai. 

model can be obtained by choosing (3i = 9/7, /32 = 7/9 
and s = 2. As we can see, again, modified gravity en- 
hances the lensing potential even if the dependence from 
Ai is different. 

III. ANALYSIS METHOD 

In order to analyze the impact of a non-standard weak 
lensing component on parameters inference we simply pa- 
rameterize the weak lensing signal by defining a fudge 
scaling parameter affecting the lensing potential power 
spectrum (see Q): 

Cf ^ ALCf. (10) 

In other words, parameter Al effectively multiplies the 
matter power lensing the CMB by a known factor. Al = 
is therefore equivalent to a theory that ignores lensing of 
the CMB, while Al ~ 1 gives the standard Icnscd theory 
and Al > 1 may suggest modified gravity or < 1. 
Since at the scales of interest the main effect of lensing 
is purely to smooth the peaks in the data, Al can also 
be seen as a fudge parameter controlling the amount of 
smoothing of the peaks (sec p6|). 

As we can see in Figure ^ the inclusion of this param- 
eter can well recover a non standard weak lensing signal. 
This is in agreement with the results presented in ||3^ 
where a one-parameter description of lensing could rea- 
sonably recover all the main information achievable from 
CMB lensed spectra. 

In particular we found that for w = —0.8 a model with 
Cs ~ 1 and Al ~ 1.06 can perfectly reproduce the same 
scenario with cj. = 0, while a modified gravity model, 
compatible with solar system bounds, can be approxi- 
mated by Al = 1.12. 

In what follows we provide constraints on Al by ana- 
lyzing simulated datasets for future experiments consid- 
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FIG. 3: Lensing potential for different dark energy (Top 
Panel) and modified gravity (Bottom Panel) models. The 
Al parameter can well approximate the higher amplitude of 
the potential and the deviations from the standard case. 



cring also B polarization modes^. Moreover, we inves- 
tigate the impact of a wrong assumption related to the 
weak lensing projected potential amplitude on the deter- 
mination of the cosmological parameters. 

We constrain the by a COSMOMC analysis of fu- 
ture mock datasets. The analysis method we adopt is 
based on the publicly available Markov Chain Monte 
Carlo package cosmomc p8|] with a convergence diagnos- 
tic done through the Gelman and Rubin statistic. We 
sample the following eight-dimensional set of cosmologi- 
cal parameters, adopting flat priors on them: the baryon 
and cold dark matter densities Wb and ujc, the ratio of 
the sound horizon to the angular diameter distance at 
decoupling, 0,, the scalar spectral index ng, the overall 
normalization of the spectrum A at fc = 0.05 Mpc^^, the 
optical depth to reionization, r. Furthermore, we con- 
sider purely adiabatic initial conditions and we impose 
spatial flatness. 



^ We ignored primordial tensor modes 
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We created several full mock datasets (temperature, E 
and B polarization modes) with noise properties consis- 
tent with Planck EPIC and §), and for SPT § 
(see Table |), assuming, as the fiducial model, the best-fit 
from the WMAP plus ACBAR analysis of Ref. Q with 
ujh = 0.0227, ujc = 0.113, ns = 0.0973, t = 0.0908 and 
Al = 1. 



Experiment Channel FWHM AT/T 



Planck 


143 


7.1' 


2.2 


faky = 1.0 








Epic 2m 


150 


5' 


0.44 


fsky = 0.85 








SPT 


150 


1' 


10 



fsky = 0.06 

TABLE I: Planck, EPIC, and SPT experimental specifica- 
tions. Channel frequency is given in GHz, FWHM in arcmin- 
utes and noise per pixel in 10"^ for the Stokes-I parameter; 
the corresponding sensitivities for the Stokes Q and U param- 
eters are related to this by a factor of %/2. 

We consider for each channel a detector noise of if^^ = 
(6*0-)^ where is the FWHM of the beam assuming Gaus- 
sian profile and a is the sensitivity IST /T both from Ta- 
ble |. We therefore add to each fiducial spectra a noise 
spectrum given by: 

ATf = w-icxp(/(Z + l)Z^), (11) 

where is given by 4 = -\/8 h\2/6. We analyze these 
datasets with a full-sky exact likelihood routine as in 
Rcf. H . In the case of SPT we also include a fjf.y ~ 0.06 
pre-factor {fsky = 1 for Planck and fsky = 0.85 for 
EPIC). 

IV. RESULTS 



Experiment 


Dataset 


Limits on Al 


Planck 


w/o Bmodes 


r.uu_Q og_g 11 


Planck 


with Bmodes 


1 p.p.+om+o.i2 

r.UU_o 06-0.12 


EPIC 


w/o Bmodes 


l.UUU_Q Qi3_Q Q26 


EPIC 


with Bmodes 


±.UUU_g Qij_Q 


SPT 


w/o Bmodes 


1 |-|f^+0.26+0.60 
i.UU_o 26-0.44 


SPT 


with Bmodes 


^^•^'-'-0.13-0.34 



TABLE II: Fiducial Model parameters and limits on Al for 
Planck and EPIC satellites and for the SPT experiment. We 
report errors at 68% and 95% confidence level. 



In Table we report the constraints on the param- 
eter achievable by future experiments. Planck will reach 



an accuracy of about ~ 6% at Icr level while EPIC will 
constrain Al at ^ 1% level (always la). It is interest- 
ing to note that the exclusion of the B-mode polariza- 
tion channel from these experiments has small effects on 
the determination of A^. In a few words, the smearing 
of temperature and polarization anisotropics induced by 
lensing provides a larger statistical evidence for Al than 
a detection of Icnscd B-mode polarization. 

This will not be the case for the SPT experiment, 
where the higher angular resolution will permit better 
measurement of the B-mode polarizations signal. In the 
SPT case, the inclusion of B modes practically halves the 
error bars on Al- However, as already discussed in Q 
the sample variance on the small scale B modes could be 
larger by a factor of ~ 10, reflecting the variance of the 
larger scale lenses that generate them. This degradation 
effect should be considered for SPT and could strongly 
affect the reported results on Al when B-modcs are con- 
sidered. Being the degradation effect much smaller for 
the TT, EE and TE modes, in what follows we only 
consider forecasts for the Planck and EPIC experiments. 

Comparing the future constraints on Al with the ex- 
pected signal from dark energy perturbations, it is clear 
that the Epic experiment may have enough sensitivity to 
constraint the sound speed parameter if w ~ —0.8. In 
order to test this, we have investigated the constraints 
achievable by Planck and EPIC on the — plane un- 
der the assumption of a fiducial model with w = —0.8, 

= and the remaining parameters fixed as above. 

In Figure ^ we plot the 68% and 95% confidence levels 
expected from Planck and EPIC with and without con- 
sidering the lensing signal. As we can see, when lensing 
is considered, is strongly constrained for larger values 
of w. Adding complementary information on w such that 
w = -0.8 ± 0.1 yields < 0.1 from EPIC at 95% c.l. 

On the other hand, the constraints obtained on Al will 
probe modified gravity scenarios. In particular, in the 
case of /(i?) models, an upper limit of Al < 1.06 from 
Planck {Al < 1.01 from EPIC) will bound the Compton 
wavelength to A? < 2 . lO^Mpc^ (Planck) and Xj < 2 ■ 
lO^Mpc^ (EPIC). Models that are consistent with solar 
system test can be significantly ruled out by EPIC. For 
Chameleons models the bound will be similar with < 
8 • lO^Mpc^ (Planck) and A? < 9 • lO^Mpc^ (EPIC). 

The accuracy reachable by Planck and EPIC will there- 
fore open the possibility of testing non standard modifi- 
cation to gravity at cosmological length scales (|Q). It 
is interesting to quantify the impact of a wrong assump- 
tion in the lensing potential on the constraints derivable 
from those experiments and the corresponding ability of 
recovering the correct solution. 

In Figures || and ^ we can see the 2-D likelihood con- 
tour plots which show the degeneracies between the lens- 
ing parameter Al and the remaining cosmological pa- 
rameters. As we can see, a degeneracy appears between 
Al and the matter density coc and the spectral index ng 
while a milder degeneracy is also present with the baryon 
density and the optical depth r. A larger value of will 
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FIG. 4: 68% and 95% c.l. constraints on the c^ — w plane from 
Planck (Red) and EPIC (Blue) when lensing is considered 
(Top Panel) and without considering lensing (Middle Panel). 
The fiducial model assumes w = —0.8 and — 0. In the 
Bottom Panel we can see the effect of adding an external 
prior on w = —0.8 ±0.05 at 68 c.l.. A combined analysis with 
CMB data could discriminate values of lower than one. 



make larger values of Ub and ns and smaller values of lOc 
and r compatible with the CMB data respect to the case 
with Al ~ 1. 

Including an uncertainty on will therefore relax the 
bound on these parameters, especially on ns and ujb- It 
is however interesting to quantify the error one could 
make if not including a marginalization over A^ when 
analyzing future datasets. 




FIG. 5: 68% and 95% c.l. 2-D Likelihood contour plots for 
Al versus different cosmological parameters from the Planck 
experiment. 




FIG. 6: 68% and 95% c.l. 2-D Likelihood contour plots for 
Al versus different cosmological parameters from the EPIC 
experiment. 



We have therefore changed the fiducial model to = 
1.3 for Planck {Al = 1.05 for EPIC) and re-analyzed 
the mock datasets either with the wrong assumption of 
Al = 1, either letting Al to vary, 

As we can see from Figures and H, not accounting for 
an higher amplitude of the lensing potential could lead 
to misleading constraints on the baryon density, the cold 
dark matter density the spectral index and the optical 
depth. In particular, if we don't marginalize over Al 
we obtain the constraints ujb = 0.02249 ± 0.00014, ojc = 
0.1150 ± 0.0013, Us = 0.9672 ± 0.0038, t = 0.0924 ± 
0.0047 for Planck and uub = 0.02264 ± 0.00006, uJc = 
0.11426±0.00034, r = 0.0931±0.0024 and Us = 0.9704± 
0.0027 for EPIC to be compared with the fiducial values 
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FIG. 7: Bias introduced in cosmic parameters inference by the 
assumption of a wrong weak lensing amplitude in the case of 
the Planck experiment. The fiducial model has Al ~ 1.3. 
The solid line show the constraints when the Ai parameter is 
let to vary, while the dotted constraints are assuming Al = 1- 
The straight black line identifies the fiducial model. 
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FIG. 8: Bias introduced in cosmic parameters inference by the 
assumption of a wrong weak lensing amplitude in the case of 
the EPIC experiment. The fiducial model has Al = 1.05. 
The solid line show the constraints when the Ai parameter is 
let to vary, while the dotted constraints are assuming Al = 1- 
The straight black line identifies the fiducial model. 



of LUb = 0.0227, uJc = 0.113, 0.973, r = 0.0908. 

When a marginalization over Al is included we recover 
the correct fiducial input parameters. 



V. CONCLUSIONS 



In this paper we have investigated the abihty of fu- 
ture measurements of CMB temperature and polarization 
angular spectra to identify a non-standard weak lensing 
signal. We have found that the amplitude of the weak 
lensing potential will be measured with a ~ 6% accuracy 
from Planck (1% from EPIC) at 68% confidence level 
This will let to test the dark energy sound speed param- 
eter c| for values of the equation of state different from 
— 1. If, for example, the dark energy of state will be con- 
strained at w = — 0.80± 0.05 at 68% c.l. then a value of 
— could be discriminated at more than 3 standard 
deviations from = 1 from an EPIC-like experiment. 

Moreover, we have found that both Planck and EPIC 
could place new constraints on modified gravity models. 
Translating the bounds on in the case of f{R) mod- 
els on the Compton wavelength Ai we have found that 
Planck wiU place a bound of Xj < 2- lO^Mpc^ while EPIC 
will reach < 2 • lO^Mpc^. For Chameleons models the 
bound will be similar with A^ < 8 • lO^'Mpc^ (Planck) 
and A^ < 9 • lO^Mpc^ (EPIC). This will present at least 
a two order of magnitude improvement on the constraint 
A^ < 4-lO'^Mpc^ from current CMB data ||9|. Finally one 
should remember that constraints on can be placed 
by considering higher order correlations and extracting 
the lensing deflection field from the CMB maps. Several 
methods have been proposed as the quadratic estimator 
method (s ee [|39| ) and the iterative method by Hirata 
& Seljak (1 40 ). Estimating the noise in the case of a 
quadratic estimator, we have found that a ~ 1% error on 
Al could be reached by EPIC of about the same order 
of the constraint presented here obtained from measure- 
ments of the angular power spectra. Since experimental 
systematics and foregrounds may strongly affect future 
CMB measurements but affect the angular spectra and 
the quadratic estimator in a different way the two meth- 
ods are complementary and will both provide valuable 
information to the dark energy component. 
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